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We explore a material design strategy to optimize the thermoelectric power factor. The approach 
is based on screening the band structure changes upon a controlled volume change. The methodology 
is applied to the binary silicides and germanides. We first confirm the effect in antifluorite Mg 2 Si 
and Mg 2 Ge where an increased power factor by alloying with Mg 2 Sn is experimentally established. 
Within a high-throughput formalism we identify six previously unreported binaries that exhibit 
an improvement in their transport properties with volume. Among these, hexagonal MoSi 2 and 
orthorhombic Ca 2 Si and Ca 2 Ge have the highest increment in zT with volume. We then perform 
super-cell calculations on special quasi-random structures to investigate the possibility of obtaining 
thermodynamically stable alloy systems which would produce the necessary volume changes. We 
find that for Ca 2 Si and Ga 2 Ge the solid solutions with the isostructural Ca 2 Sn readily forms even 
at low temperatures. 

PACS numbers: 


Despite their importance, the discovery of new mate¬ 
rials are often based on trial a nd e rror. High-throughput 
(HT) computational screenin^^^ is an important step 
towards identifying materials with desired properties in 
a more systematic way. Thermoelectric (TE) materials 
are attractive for such computational searches because 
continuous development of computational methodology 
means that all parts of the TE figure of merib zT can in 
principle be calculated from first principlesIn prac¬ 
tice computational HT searches for new TE materials 
have focused on parts of the zT = S’^aT/n, where S is 
the Seebeck coefficient, cr the electrical and k the ther¬ 
mal conductivity.^^^ Despite this, there are now a few 
works where computational screening has led to high 
performance TE materials that could be experimentally 

realized. 

Beyond screening known compounds, there still ex¬ 
ists great challenges designing new materials with spe¬ 
cific properties. This is especially the case for electronic 
structure dependent properties, which have highly non- 



FIG. 1: Schematic illustration of volumetric band alignment 
for a n-type material. VBM stands for the Valence Band 
Minimum. 


trivial dependencies on the atomic structure.!^ Eor TE 
materials one strategy for designing new alloys with op¬ 
timized properties is by a controlled volume change. We 
label this procedure volumetric band-structure alignment 
(VBA). The idea is illustrated in Fig.[^ where the energy 
dependence of two bands vary differently upon a change 
of volume. Thereby a scenario can occur when the band 
edges are aligned, as schematically illustrated in the mid¬ 
panel of Fig. [l] How this optimizes the TE power factor, 
PF = S'^a, can be understood by considering the gener¬ 
alized transport coefficients, 

£(“) = ^2^ cr(e)(e_^)“^_^^d£, (1) 


where / is the Fermi-distribution and (T(e) the transport 
distribution. The Seebeck and electric conductivity are 
given as S = /qT and cr = respectively. For 

two channels (labelled ' and ") conducting in parallel, 

£(a) 

is given as the sum of the contributions from each 
channel, so that the PF is. 


S'V = 


1 (£(!)'+£(!)")" 
q2rp2 £( 0 )' £( 0 )" 


( 2 ) 


Without loss of generality we can write 
and . Thereby it is clear that the PF will 

be increased if, 


(1 -|- Oi)^ > 1 -|- ctQ. (3) 

and thereby S will be significant when the chemical 
potential is located within the relatively narrow energy 
window defined by df/de of a band edge. Aligning band 
edges with similar effective masses within df/de, will 
satisfy Eq. (§ and lead to an increased PF when com¬ 
pared to the largest of the individual contributions. Such 
bandstructure characteristics are found in several mate¬ 
rials with complex carrier pocket shapes which results in 
enhanced TE properties 
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VBA has been successfully applied to optimize the PF 
in the n—doped Mg 2 Sii_xSna; ( and M g 2 Gei- 3 ;Sna;) and 
in p—doped PbTei_ 2 ;Sea; alloys.l23^ In both the cases 
the volume was controlled by alloying. Consider as an 
example the alloy composition MXi_ 2 ;X’a;. According to 
Vegard’s law, the volume of the resulting alloy system, 
will be directly proportional to the fraction of the new 
composition added, i.e. x. The volume will also be influ¬ 
enced by the thermal expansion, and this together with 
alloying offers the possibility to control the peak temper¬ 
ature of TE performance.!^ 

The question whether it is possible to optimize the 
band structure by VBA for a given compound, is very 
difficult to answer from intuition alone. It depends on 
the detailed band structure and how the different bands 
react to a volume change. The effect of VBA can be 
strongly affected by changes to the band gap and doping 
level. Furthermore, it can be very difficult to predict 
whether a given compound can be alloyed or not. E.g. 
while a solid solution of Si and Sn in Mg 2 Sii_a;Sn 3 ; can 
form even at low temperature^^lHlZ] it jg -y^rell known that 
Sn is hardly soluble in diamond-Si. 

The idea behind the present work is to explore a HT 
computational strategy to identify systems where VBA 
can be applied. The approach is based on screening the 
volumetric effect on the band structure and calculate ab- 
initio thermodynamics to asses the possibility of alloying. 
It is shown that systems allowing VBA are quite rare and 
that the procedure leads to a very strong screening of po¬ 
tential candidates. We focus on the electronic part of the 
zT. However, alloying will also be advantageous in terms 
of reducing the thermal conductivity either by lowering 
mass-disorder scattering or by usin g n atural solubility 
limits to nano-structure the material.!^ 

We will explore binary alloys of the group 14 ele¬ 
ments. Silicon is the second most abundant element on 
the earth’s cruslP^ making such systems quite attractive. 
Stable alloys and solid solutions with optimized electronic 
structures and transport properties have been realized 
with many binary silicides plus the corresponding ger- 
manides or stannides compounds.!^ While Ge and Sn 
are not as abundant and cheap as Si, they are certainly 
not rare elements either.!^ Other crucial advantages that 
these M-Si/Ge/Sn alloy systems enjoy, is the ease with 
which they can be doped. Indeed silicide group com¬ 
pounds exhibit a high density of mobile charge carriers 
upon doping of up to n = 10^° — 10^^ cm“^. Further¬ 
more, they often have high melting points making them 
attractive stable candidates for high-T TE.I^ 


I. METHODS 

A. High-throughput scheme 

Our HT scheme is illustrated in Fig. We consider 
binary M-X systems, where M is a metal taken to be a 
3d, id or 5d Transition Metals (TM) or 2s, 3s and 4s 
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FIG. 2: Flowchart of the HT procedure employed to inves¬ 
tigate the candidates for volumetric band engineering. The 
criterion for good thermoelectrics, i.e. zT > 0.4, within our 
scheme is checked at T = 600 K and n = 2 x 10^° cm“® with 
zT evaluated from Eq. 0- 


Alkali (A) or Alkaline Earth (AE) metals and X is Si, Ge 
or Sn. The first step of our work is to identify the stable 
structures in the constituent binary system. The phase 
stability is evaluated within a l^h throughput formalism 
based on our previous works.!^ For each metal silicide 
(or germanide) combination, the crystal structures were 
relaxed and the formation energy was calculated. The 
candidate structures were generated by the extracting 
the known structures from the Pearson’s databas^^ and 
appropriate substitution of the known structures with 
similar atoms (in same group of Periodic Table). Further¬ 
more, we have also incorporated the silicides structures 
from our previous worlP^ and consequent substitution of 
similar atoms. 

In total, we have investigated 3150 different com¬ 
pounds. For a given M-X binary system, we select the 
structures that have the difference in formation energies 
from the corresponding convex hull, AEh, lower than 
50 meV/atom. The AEh < 50 meV/atom tolerance has 
statistically been shown to contain 80 % of the experi¬ 
mentally known compounds in the TM-Si systerrP^ and 
narrows the original number to 522 thermodynamically 
stable compounds, Fig.|^ 

Thereafter, the self-consistent calculations for these se¬ 
lected structures were performed using the (L)APW-l-lo 
methocP^ implemented within the WIEN2k code.!^. 
These were followed by bandstructure calculations on a 
finer k-mesh of k-points in the full BZ. All cal- 

culations in this work are reported for Perdew-Burke- 
Ernzerhof (PBEjP^ exchange-correlation potential. Sub¬ 
sequently, the electronic transport properties were eval¬ 
uated using the BoltzTrap code.!^ The code evaluates 
the Seebeck coefficient on an absolute scale and electri¬ 
cal conductivity in terms of the carrier relaxation time 
r using the rigid band approximation. To evaluate the 
potential for TE energy conversion we utilize a modified 
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version of the definition of figure of merit, zT, 


ing Gibbs free energy is, 


zT = 


S2^T 


( 4 ) 


RT 


H-[a;lna; + (1 — cc) ln(l — a:)l (6) 

n + 1 '■ ■' 


The denominator in Eq. 0 accounts for the total ther¬ 
mal conductivity of the material, with its electronic part, 
Kei = Lo^T, written according to the Wiedemann-Franz 
law. Here Lq is the Lorenz number, Lq = 2.44 x 10“® 
WHK“^. Finally the phononic part of thermal conduc¬ 
tivity is Kph. All the transport quantities in Eq. Q are 
extracted on an absolute scale, apart from ^ which is 
set at 10^"* W ms“^.I^In reality, the magnitude of 
will be dictated by both intrinsic properties such as the 
phonon band structure and phonon- and electron-phonon 
coupling and by the details of sample preparation control 
e.g. grain size (nanostructuring) or the existence of mul¬ 
tiple phases. zT, as defined in Eq. Q, epitomizes the 
electronic contribution to the TE performance a nd has 
been validated in the prediction of TE materials 
It can be viewed as a descriptor for identifying poten¬ 
tial high TE performance under the assumption of a low 
thermal conductivity, a long electron life time and the 
possibility of reaching the optimal doping. Setting a cri¬ 
teria of zT > 0.4 [evaluated using Eq. 0] at T = 600 K 
and 2 x 10^° cm“^ reduces the number of potential can¬ 
didates to 29. 

The procedure was then repeated at unit-cell volumes 
in the ±7% range of the equilibrium volume. The volume 
dependence of the descriptor, i.e. zT{V), is evaluated 
and we identify zT^pt and zTq which are the magnitude 
of zT at the volume which maximizes zT(V) and at the 
calculated equilibrium volume respectively. Thus, for any 
given material, if > 1-1, it will be of interest for 

the purpose of this work. As can be seen in Fig. [^this 
lowers the number of candidate structures to only eight. 

The Hnal step is to consider whether the volume change 
by alloying is thermodynamically feasible. Consider a 
mixture of 1 — a; mole fractions of the binary ABn and 
X mole fractions of the binary ACn producing the alloy 
ABn-xCx- The excess energy which is required to obtain 
the alloy, also referred to as the mixing energy, is 


Aii/niix(a;; AB„, AC„) = AEf{ABn-xCx) 

xAEf{ACn) + {n — x)AEf{ABn) 
n 

( 5 ) 

where AEf are the formation energies per atom for the 
given compounds. We have calculated the mixing ener¬ 
gies from Eq. ([^, by taking Special Quasi-random Struc¬ 
ture (SQS) alloy distributions for different values of x us¬ 
ing the ATAT codJ^. For solids, since at ambient pres¬ 
sures, PAxaixV term is negligible, the enthalpy of mixing 
can be approximated as AiLmix = AE^nix- Thus the mix- 


where, the last term accounts for entropy of mixing. 
For each of the alloys, the minimum temperatures, at 
which the entropy gain in alloy formation compensates 
the energy cost of mixing, can be obtained by minimiz¬ 
ing Eq. Q. We will report T^, which is the maximum 
temperature of the boundary of corresponding miscibil¬ 
ity. As will be discussed later, this in the end limits the 
number of potential candidates to only four. 


B. Identifying Mg2X as a promising candidate. 

In the following we will introduce the HT scheme 
by example, using the Mg 2 Si/Ge/Sn system, which is 
well known for its potential VBA.I^ In the top panel 
of Fig. 1^ the temperature dependance of the enhance¬ 
ment of the TE properties upon volume optimization, 
i.e. PE{Vopt) — PE{Vo) and — «S'(Vb) for the anti¬ 

fluorite structure of Mg 2 Si (Fig.[^) and Mg 2 Ge (Fig. Ep) 
are shown at a doping of n = 2 x 10^° cm In the 
bottom panel Fig. 0, the bandstructures for Mg 2 Si at 
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FIG. 3: Volumetric Band Alignment in Mg 2 X. In the top 
panel. Fig. a and Fig. b respectively show the change in 
BE(Vopt) - PF{Vo) (in red circles) and ^(Vopt) - S{Vo) (blue 
triangles), as a function of temperature at n = 2 x 10^° cm“®. 
Furthermore, in both the graphs the dashed lines correspond 
to n—doping while the bold line represents the data for the 
p—doping scenario. The bottom panels. Fig. c and Fig. d re¬ 
spectively show the bandstructures for Mg 2 Si and Mg 2 Ge, at 
Vo (red) and Vopt (black). The insets represent the zoomed-in 
view of the conduction band minimum. The conventions used 
in this figure will be used in the forthcoming figures presenting 
similar results for other candidates in this paper. 


S(V.„)-S(V„) [,.V/K] 
































4 



v(A^) 

FIG. 4: Formation energy vs volume for Mg 2 Xi_a;Yx alloys. 
Here X,Y are Si, Ge and Sn. The insert shows the miscibility 
gap calculated by minimizing Eq. §. 


equillibrium volume and the optimized volume are plot¬ 
ted. Likewise, Fig. w illustrates the bandstructures at 
Vo and Vopt for Mg 2 Ge. 

In Fig. 1^ and Fig. we observe that for the p—type 
behavior the change in PF with volume is negligible for all 
T. On the other hand, for n-type behavior, the PF shows 
a pronounced increase at the optimized volume for both 
the compounds, that further enhances with T. Moreover 
S also exhibits an increase. Note that for n—doping, the 
higher the negative value of S{Vopt) — S{Vo), the more 
is the desired enhancement of the thermopower and zT. 
Also, it is well known that Mg 2 Si exhibits a persistent 
n—type conductivity under any kind of crystal growth 
conditionsWhen we analyze the variation in the 
conduction bands (CB) of the compounds with volume 
(Fig. and Fig. |^), the large increment in the PFs for 
n—doped Mg 2 Si and Mg 2 Ge becomes quite apparent. At 
the X point in the BZ near the conduction band mini¬ 
mum (GBM), the first and second GBs directly coincide 
at Vopt for Mg 2 Si and are only a few meV appart for 
Mg 2 Ge (see insets). Therefore, even at low dopings, the 
charge carriers (electrons) residing in both the bands will 
contribute to enhance the PF and Seebeck coefficient, as 
opposed to the contributions coming from electrons in 
only a single band, at equillibrium volume. 

For Mg 2 Si, at a doping of n = 2 x 10^° cm ^ and 
T = 600 K, the highest increment in zT is observed at 
volume increase of 2% from the equilibrium volume, Ta- 
ble|Tj According to Vegard’s law this corresponds roughly 
to alloying Mg 2 Si with 10% Mg 2 Sn. This value will also 
be dependent on the carrier concentration and temper¬ 
ature at which the transport properties are calculated. 
The important rule is if a large VBA effect is observed at 
a small volume change that in principle can be attained 
by alloying. The second crucial criterion is whether the 
alloying is thermodynamically feasible. The energy of 
formation vs volume for Mg 2 Xi_ 2 ;Y 2 ; solid solution is de¬ 


Gompound 

zTopt 

zTopt/zTo 

AEh [AEt] 
(meV/atom) 

doping 

^opt 

(%) 

MgzSi 

1.07 

1.33 

0[0] 

n 

2.0 

Mg2Ge 

1.17 

2.19 

0[0] 

n 

5.0 

Ga2Si 

0.80 

1.33 

0[0] 

n 

5.0 

Ga2Ge 

0.81 

1.25 

0[0] 

n 

2.5 

GagGes 

0.96 

1.47 

37.4[17.9] 

p 

6.1 

d-MoSi2 

0.77 

1.75 

27.3(180.8] 

n 

3.0 

o-Fe2Ge3 

0.47 

1.19 

0.1(22.5] 

n 

3.0 

t-Fe2Ge3 

0.54 

1.38 

0]22.8] 

n 

3.0 


TABLE I: Summary of eight candidates that exhibit the VBA 
effect. We summarize zTopt and zTopt/zTo magnitudes at 
doping of n = 2 X 10^° cm“® and T = 600 K. The distance 
from the convex hulls, AiJh [AAf"] of the candidate [alloying 
possibility], carrier sign and Vopt (expressed as an increment 
from Vq) are also summarized. 

picted in Fig. Here X is Si or Ge while Y is Ge or Sn. 
It can be seen that the excess mixing energy required for 
the formation of stable alloys is small in magnitude. As 
a result the optimal Sn content can be reached at mod¬ 
erate growth temperature, Fig. |^nsert. AE^^ix and 
are in good agreement with previous DFT calculatiorl^, 
especially considering that we employ a somewhat differ¬ 
ent computational approach using SQS based supercells 
and linear interpolation between the calculated points 
when minimizing Eq. ([^. We also find a good agree¬ 
ment with GALPHAD and experimental results for this 
system .1231^. 

Thus in the particular cases of antifluorite Mg 2 Si and 
Mg 2 Ge, the VBA can be conveniently achieved through 
alloying. While Mg 2 Sii_a;Gea; alloy thermodynamically 
exists as a solid solution even at low temperatures, for 
Mg 2 Sii_ 2 ,Sn 2 , and Mg 2 Gei_a;Sn 2 , alloys the desired vol¬ 
ume increment can be achieved at the expense of only a 
small magnitude of mixing energj®ll2l 


II. RESULTS AND DISCUSSION 
A. Electronic screening 

Based on the descriptors discussed above other can¬ 
didates for VBA effect will now be identified. The top 
and bottom panels in Fig. illustrates xTopt vs zTq re¬ 
spectively for all silicides and germanides, at a doping of 
n = 2 X 10^° cm“^ and T = 600 K. Each data point in 
Fig. [5] represents a particular compound. 

Thus Fig. [^contains three vital descriptors to identify 
candidates for VBA, (i), the magnitude of zTq, which 
indicates if a particular candidate is a good TE mate¬ 
rial, (ii), the value which helps us to identify the 

promising candidates exhibiting the volumetric enhance¬ 
ment of their TE properties, (iii) the magnitude of the 
distance from the corresponding convex hull of the possi- 
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FIG. 5: zTopt vs zTo for all silicides (top) and germanides 
(bottom) at doping n = 2 x 10^° cm“® and T = 600 K. Fur¬ 
thermore, the scatter points are colored based on the dis¬ 
tance from the convex hull, AFif", of the corresponding Sn 
compound with the same structure as the host binary. For 
example, in the case of the CagGes compound, Ai?h of the 
corresponding iso-structural GagSns is used for the color code, 
and so on. The candidates that exhibit VBA, > 1.1, or 

have a large value of zTq are labelled. Also, if a compound 
crystallizes in more than one phase, the corresponding struc¬ 
ture is also indicated with a prefix o- for orthorhombic, t- for 
tetragonal and af- for anti-fluorite phase. 


ble alloying choices, which would produce the necessary 
volumetric change upon alloying (according to Vegard’s 
law). 

A summary of the promising silicide and germanide 
structures that exhibit a large VBA effect are shown in 
Table 1^ with the corresponding value of zTopt, zT^^tl zTq 
and Ailh [Ai?®"]. Eight such compounds were identified 
to exhibit encouraging VBA effects upon using the mag¬ 
nitude of zTopt/zId and AE^ as descriptors. 

In general, with increasing doping we observe that 


Alloys 

^alloy 

AFmix(a: = 0.25) 
(kj/mol.atom) 

(K) 

X 

(T = 800 K) 

MggSii-xSux 

0.09 

1.997 

1281 

0.171 

MggGei-xSux 

0.26 

1.595 

1023 

0.138 

C3/2 Sil — xSllx 

0.30 

0.013 

78 

all 

CagGei-xSux 

0.18 

0.015 

97 

all 

CSjQ G 65 — xSllx 

0.28 

4.495 

> 5000 

0.019 

/3—MoSig-xSrix 

0.07 

30.980 

> 5000 

0.008 

o-FegGea-xSux 

0.10 

32.746 

> 5000 

0.004 

t-FegGea-xSux 

0.10 

29.681 

> 5000 

0.005 


TABLE II: Alloy thermodynamics of the candidates. Xaiioy is 
the amount of alloying component, calculated using Vegard’s 
law, required to obtain the volume change for the VBA effect. 
AEinix(x = 0.25) is the mixing energy required to form the 
alloy (at x = 0.25). Tm is the maximum temperature of the 
boundary of the miscibility gap and x(T = 800 K) is the 
amount of alloying component which can be accomodated at 
T = 800 K in the alloy. 


more number of candidates show the VBA effect. This 
is due to a large value of zT itself that is improved upon 
doping. However the best candidates, as listed in Ta¬ 
ble exhibiting VBA remained unchanged. Apart from 
Mg 2 Si and Mg 2 Ge, discussed previously, orthorhombic 
Ca 2 Si and Ca 2 Ge, hexagonal CagGes and hexagonal- 
MoSi 2 show encouraging results. 

Note that in Fig. there are compounds that show 
very little to no VBA effect but exhibit large zTq. In¬ 
deed, despite our semi-emperical determination of zT, we 
correctly predict most of these binaries with large values 
of zTq, that are already established as encouraging TE 
materials.!^ Compounds which exhibit zTg in excess of 
1.0, at n = 2 X 10^° carriers per cm“^, are the orthorhom¬ 
bic and tetragonal structures of Ru 2 Ge 3 , cubic RuSi (at 
low temperatures) and tetragonal NaGe and NaSi. Fur¬ 
thermore, the following compounds have zTg between 0.5 
and 1.0 (not listed in Fig. [^: hexagonal CrSi 2 , Mn 4 Si 7 , 
FeSi 2 , ReSii. 75 , and CoSi. 


B. Feasibility in alloying 

We will now explore the posibility of forming thermo¬ 
dynamically stable alloys among the candidates for VBA, 
identified based on the descriptors discussed so far. These 
alloys are summarized in Table. |TT] together with pro¬ 
portion of alloying component Xaiioy required to produce 
the highest optimized PFs in the candidates. Note that 
the optimal VBA for most candidates are achieved at 
X < 0.30. 

In the previous section we have established that Ca 2 Si 
and Ca 2 Ge are attractive candidates for VBA effect 
based on the descriptors which characterize their trans¬ 
port properties. In Fig. we illustrate the energy of 
formation of orthorhombic G& 2 ^i-x^x solid solutions as 
a function of volume. We observe that for x = 0.25 
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FIG. 6 : Formation energy vs volume for Ca 2 Xi_ 2 ;Ya; solid 
solutions (X and Y are Si, Ge and Sn). 




C 



FIG. 7: Ordered alloy structure is shown here for a 2 x 2 x 2 
supercell. Here the blue balls are Ca atoms, the grey balls 
are X atoms and the magenta ones are Y atoms. Note that 
an unitcell contains 12 atoms, i.e. Ca 8 X 2 Y 2 composition. 
The shortest crystallographic axis, b, is perpendicular to the 
plane. This structure is related to the Go 2 Si (or TiNiSi) struc¬ 
ture type with space group PnmcP^. 


and X = 0.5, the SQS generated random structures are 
actually lower in energy in comparision to the parent 
structures. We therefore generated all the possible or¬ 
dered structures of Ca 64 Si 24 Sn 8 and Ca64Sii6Sni6 within 
a 2 X 2 X 2 supercell. We found a previously unreported 
ordered structure in which the Sn atoms order along the 
short b-axis to be marginally lower in energy in com¬ 
parision to the random SQS structure chosen. Fig. [ 7 ] 
illustrates this ordered structure, which can be thought 
of as being stabilized by Sn-Sn bonding along the short 
axis of the orthorhombic cell and by the shortest Sn- 
Sn distance within the unitcell. In Table. |ll] the energy 
difference between the ordered and disordered alloy is re¬ 
ported as APmix. The energy difference is very low and 
the alloying necessary to obtain the desired volume in¬ 


crement should be readily achieved through Ca 2 Sii_xSnx 
and Ca 2 Gei_xSnx alloys. Thus Ca 2 Si and Ca 2 Ge com¬ 
pounds are confirmed to exhibit the VBA effect and the 
required volume change can also be conviniently achieved 
through alloyi:^ Interestingly, Ca 2 Sn has been previ¬ 
ously proposecP^l as a thermoelectric material. 

Unfortunately for all the other alloys, we obtained the 
mixing energy AUmix to be very large (see Table. 0 - 
Therefore for GagGes, /S-MoSig, o-Fe 2 Ge 3 and t-FegGea, 
the optimal VBA volume cannot be expected to be at¬ 
tained by the simple alloy choices listed in Table. |T^ 

Interestingly, we observe that for all four alloys having 
a large AAmix(a; = 0.25) in Table |nj the correspond¬ 
ing pure compounds are found above the convex hull. 
Table |Ij E.g. in the case of /3-MoSi2 we had obtained 
that /3-MoSn2 is far from the convex hull (i.e. « 180 
meV/atom). This could indicate that AA®" in Fig. 
can serve as a minimal requirement for the possibility of 
forming a stable alloy system. Such a descriptor could 
help to avoid the computationally more expensive SQS 
supercell calculations. 


C. Origin of high zT(Uopt) in enconraging structures 

In the following, we shall investigate the role played by 
the electronic structure on VBA effect for the candidates. 
The main focus will be put on GaSi 2 and CaGe 2 where 
it has been shown that the necesarry volume expansion 
should be achieved by alloying. We will furthermore dis¬ 
cuss /3-MoSi2 and GagGeg. For these candidates, while 
the optimal VBA volume can not be attained by alloying, 
we point out that even a small volume increase generated 
from either alloying or thermal expansion could improve 
their transport properties through the VBA effect. 


1. Orthorhombic Ca 2 Si and Ca 2 Ge 

Both CagX compounds show a large volumetric en¬ 
hancement in their TE properties. The top panel in Fig.[^ 
illustrate the temperature dependance of the difference in 
PF (and S) between I4pt and Vg, for GagSi in Fig.and 
CagGe in Fig.|§). Likewise, the bottom panel shows the 
bandstructures of the two compounds (Fig.[^ for GagSi 
and Fig. w for GagGe), at the two volumes. 

The p—type behavior for both the compounds show a 
negligible amount of increase in PF and S at I4pt. Ac¬ 
cordingly, the VBM, around the vicinity of F point, is 
relatively unchanged. The situation is different for the 
CBM which has contributions coming from the F and T 
points at Vq- Due to this we have = 0.60 and 

^yCasGe ^ q at T = 600 K and 71 = 2 X lO^o cm-^. 
At optimal volume, n-type Gag Si and CagGe show an 
increase in their TE properties to = 0.80 and 

= 0.81. From the bandstructures of CagSi in 
Fig. [^, one can understand this volumetric enhance¬ 
ment. There is a lowering of the CBM at the T and 
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a) b) 





FIG. 8: The temperature dependance of PF{Vopt) — PF{Vo) 
and 5(I4pt) —S'(Vo) for Ca 2 Si and Ca 2 Ge respectively in Fig. a 
and Fig b. Furthermore, in Fig. c and Fig. d we depict the 
corresponding bandstructures (Vopt in black and 14 in red; 
Fig. c for Ca 2 Si and Fig. d for Ga 2 Ge). The orthorhombic 
structure for these compounds are thermodynamically stable. 
Both the n—doped compounds exhibit an increase in PF and 
S that can be justified by the increased DOS around the CB 
manifold. Note also the close band alignment of the first two 
GBs at the F point in BZ, for Ga 2 Ge. 



FIG. 9: The BZ of orthorhombic Ga 2 Si at I4pt, illustrating 
the electron pockets (in blue) at the high symmetry T and 
F points as well as the low symmetry pocket around Y. The 
electron pocket around the Y point does not exist at Vq. 


around the Y point. This aligns several pockets close to 
the CBM edge, which is further illustrated in Fig.|^ The 
constant energy surface of the lowest conduction band of 
Ca 2 Si at Vopt shows the electron pockets at the T and T 
points as well as those along the T — Y lines. The lat¬ 
ter electron pocket does not exist at Vq. In the case of 
Ca 2 Ge, a similar lowering of the lowest conduction band 
along the T — Y lines is observed in Fig. at Vjpt- Sim¬ 
ilar carrier pocket shapes have also been identified as the 
source of high a in isostructural orthorhombic Ca 2 Pb and 

Sr2Pb.ll^ 


2. Hexagonal CagGes 


a) 



Ll_ -'-'-'-'-'-'- ^ 

CL 250 300 350 400 450 500 550 600 


T(K) 



FIG. 10: T depeirdence of PF{Vapt) - PF{Vo) and S(Yopt) - 
S{Vo) for the hexagonal GagGes (top). The corresponding 
bandstructures at the two volumes are shown in bottom panel. 
For CagGes, there is an opening of the bandgap with volume 
change that is responsible for the increase in its PF and S. 


We shall now discuss the hexagonal structure of 
CagGes which also exhibits an increase in its TE prop¬ 
erties with volume change. In the same fashion as the 
data presented in the previous sections. Fig. [TO^ shows 
the temperature depedance of the change in PF and S at 
the two volumes, Vjpt and Vq. The bottom panel shows 
the bandstructure at the two volumes. 

The situation for CagGes is quite interesting. Upon an 
increase in volume, there is an opening of the bandgap 
at the F point in its BZ. This metal to semiconduc¬ 
tor transition with volume change is the sole reason be¬ 
hind the enhancement of PF and S at V)pt, graphed in 
Fig. 10 D. Interestingly, a very high = 0.96 and 

zTop^JzTo = 1.47 (at T = 600 K and n = 2 x 10^° cm“^) 
was observed. 


S. Hexagonal MoSi 2 

Molybdenum di-silicide, MoSig, exists as a—MoSi 2 , 
having a tetragonal body-centered packing with a space 
group of lA/mmm and as /?—MoSig that has a hexago¬ 
nal closed packing arrangement with space group P 6222 . 
Both the structures are composed of Mo and Sig lay¬ 
ers in which Mo atoms are surrounded by 6 Si atoms. 
Both the allotropes of MoSig exhibit properties such 
as high melting points, low resistivity and high me¬ 
chanical strength. While /3—MoSig is found to be 27.3 
meV/atom above the convex hull, it can be conviniently 
synthesized from M o and Si powders using spark plasma 
sintering techniqued^^^^, thus confirming the AP/i < 


S(V„p,)-S(V„) [,xV/K] 

























































T(K) 



FIG. 11: We depict the temperature dependence of 
PF{Vopt) - PF{Vo) and S(l/opt) - S{Vo) for the /S-MoSh 
in Fig. a. In the bottom panel (Fig. b), we present the band- 
structures at the two volumes {Vopt in black and 14 in red) 
for /3—MoSi 2 . n—type MoSi 2 shows an appreciable increase 
in both PF and S at 14pt upto T = 350 K. This increase 
can be explained by the lowering and aligning of the first CB 
manifold along multiple directions in its hexagonal BZ. 



FIG. 12: The BZ of (?-MoSi 2 illustrating the Fermi surface 
at Fopt. We observe electronic contributions connecting sev¬ 
eral low symmetry points accross the BZ in the Fermi furface 
which is responsible for the VBA effect (also see Fig. 


50 meV/atom criteria. Moreover, since a—MoSi 2 , which 
is extensively used in microelectronics, can be both p and 
n doped effectively, the same could be true for the less 
studied (3 phase. 

Fig.[TT}i shows the T dependance of difference in the 
thermoelectric properties {PF and S), between t4pt and 
Vq for /3—MoSi 2 . Additionally, Fig. [njr shows the band- 
structures at the two volumes. For n-doped MoSi 2 both 
the PF and S increases drastically by a volume expan¬ 
sion up to T = 350 K. The enhancement of PF and S 


is steadily reduced at temperatures above 350 K due to 
the small band gap. Please note that the results for all 
the candidates in Table. |T] are presented at T = 600 K. 

The explanation for the observed magnification in PF 
and S with volume, can be pinned down to the lowering in 
energy of the first CB, along multiple directions in its BZ 
observed in Fig. [n)3. These changes in the bandstructure 
not only causes an increased DOS of carriers around the 
CBM, but also an increased number of electron pockets. 
The Fermi surface demonstrating the electronic contri¬ 


butions at Fopt under doped scenario is shown in Fig. 12 
Both the above discussed factors improve PF and S'for 
n-doped MoSi 2 with volume enlargement. Consequently 
at a doping of n = 2 x 10^° cm“^, we obtain the observe 
the largest VBA effect in this work of = 1.07 

and zTopt/zTo = 3.84 (at T = 450 K) for n doping sce¬ 
nario. The same values at T = 600 K are tabulated in 
Table. U 


III. CONCLUSION 

In this paper a computational HT-scheme to identify 
compounds where the thermoelectric properties can be 
optimized by alloying is presented. 

We confirm that Mg 2 Si and Mg 2 Ge exhibit large en¬ 
hancement of their thermoelectric properties with vol¬ 
ume. We report for the first time that Ca 2 Si and Ca 2 Ge, 
hexagonal MoSi 2 and CagGea could exhibit increased 
thermoelectric properties due to a volumetric band align¬ 
ment. In the cases of Ga 2 Si and Ca 2 Ge the solid-solutions 
with Sn can be obtained at the expense of a negligi¬ 
ble amount of mixing energy and therefore the volume 
changes can be obtained by alloying. Among the re¬ 
maining candidates, the volume increase for the VBA 
effect is thermodynamically difficult to achieve by alloy¬ 
ing. However, we have obtained that /?—MoSi 2 exhibits 
a significant increase in its thermoelectric properties due 
VBA, by the virtue of favorable changes to its bandstruc¬ 
ture with volume. Thus, in this case even a small vol¬ 
ume change due to alloying or thermal expansion may 
enhance its thermoelectric properties. Finally, we have 
established that Ai?®" can be a reliable descriptor to 
provide initial information on the stability of the corre¬ 
sponding alloys. 

We have focused on systems where the electronic PF 
can be optimized with controlled volume changes. Al¬ 
loying will also decrease the lattice part of the thermal 
conductivity, which will further improve the thermoelec¬ 
tric performance of the candidates. 
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